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Q \ Cluster decays of rare-earth nuclei are studied with a view to look for neutron magic shells 

■ for the 5oSn nucleus as the daughter product always. The ^"^"Sn and ^^^Sn radioactivities 

, P^ \ are studied to find the most probable cluster decays and the possibility, if any, of new 

neutron shells. For a wide range of parent nuclei considered here (from Ba to Pt) ^^C from 

^ \ ^^^Ba and Ni from ^^"^Pt parent are predicted to be the most probable clusters (minimum 

decay half-life) referring to ^*^°Sn and ^^^Sn daughters, respectively. Also, ^^Mg decay of 
^^^Sm is indicated at the second best possibilty for ^'^"Sn-daughter decay. In addition to 
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these well known magic shells (Z=50, N=50 and 82), a new magic shell at Z=50, N=66 
(^^®Sn daughter) is indicated for the ^^Ni decay from ^^°Pt parent. 
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1 Introduction 

Since the discovery of ^^C-decay from ^^^Ra by Rose and Jones [1] in 1984, many other ^^C- 
decays from other radioactive nuclei (^^^Fr, 221,222,224,226 j^g^^ 223,225 ^^ ^^^ '^'^^Th) and some 
12 to 13 neutron-rich clusters, such as ^°0, ^^F, ^^'^^"^^Ne, ^s.soMg and ^^'^''Si, have been 
observed experimentally for the ground-state decays of translead ^^^Th to ^^^Cm parents 
[2-5], which all decay with the doubly closed shell daughter ^'^^Pb (Z=82, N=126) or its 
neighboring nuclei. Theoretically, such an exotic natural radioactivity of emitting particles 
(nuclei) heavier than a-particle was already predicted in 1980 by Sandulescu, Poenaru 
and Greiner [6] on the basis of the quantum mechanical fragmentation theory (QMFT) 
proposed by [7,8]. Todate, '^^Si is the heaviest cluster observed with the longest decay 
half-life ever measured {logiQTi/2{s) = 29.04) from ^ss-jj p^j^gj^t [9]. Recently, Poenaru 
et.al., extended the region of possible emitted clusters Ac = 14 — 34 measured in the 
region of emitters with Z = 87 - 96 to superheavy elements up to 124 [10]. In this 
systematic heavy particle radioactivity they consider not only the emitted cluters with 
atomic numbers 2 < Zc < 29 but also heavier ones up to Zc = Z — 82, around ^"^^Pb a 
doubly magic daughter. For this purpose they used Analytical Superasymetric Fission 
(ASAF) model and estimated the half- life for ^^^Sn emission from ^^^Fm (Q-value = 
252.129 MeV) and for ^^oxe emission from ^^^Ri (Q-value = 274.926 MeV): logwT^'^^s) 
= 4.88 and IoqiqT^^ {s)= 0.53, respectively. They are in agreement with experimental 
values for spontaneous fission [11]: 4.02 and 0.32, respectively. 

Keeping in mind the doubly magic nature of the ^°^Pb daughter, a second island 
of heavy-cluster radioactivity was predicted on the basis of analytical superasymmet- 
ric fission model (ASAFM)[12] and preformed cluster model (PCM) [13], in the decays 
of some neutron-deficient rare-earth nuclei in to ^"^"^Sn (Z=N=50) daughter or a neigh- 
boring nucleus. Furthermore, Kumar et al., [13] predicted another doubly closed ^^^Sn 
(Z=50,N=82) daughter radioactivity, for decays of some selective neutron-rich rare-earth 
nuclei. More ecently, an unexpected increase in E2 strengths has been reported between 
the midshell isotope ^^^Sn (Z=50, N=66) and its lighter neighbor, ^^^Sn [14], also a new 
shell closure at N=90 is predicted for the ^^"Sn isotope on the basis of shell model cal- 
culations [15]. Experimentally, several unsuccessful attempts [16-19] have been made 
to measure the ^°°Sn-daughter radioactivity from the ^^^Ba parent nucleus produced in 
^*Ni-|-^*Ni reaction. Instead, a new phenomenon of intermediate mass fragments (IMFs, 



with 3 < Z < 9), also referred to as 'clusters' or 'complex fragments', emitted from the 
excited compound nucleus, was also observed [20] . It is worth mentioning that intermedi- 
ate mass fragments are mostly found in reactions at intermediate incident energies where 
colliding nuclei breaks into many pieces[21]. 

In this paper,heavy cluster emissions of rare-earth parents (329 cases) with soSn always 
as the daughter product is considered. The new experimental mass table [22] and the 
theoretical masses [23] are used to determine the released energy. Specifically, emission 
of various isotopes of C, O, Ne, Mg, Si, S, Ar, Ca, Ti, Cr, Fe and Ni, are considered 
respectively, from neutron-deficient to neutron-rich Ba, Ce, Nd, Sm, Gd, Dy, Er, Yb, Hf, 
W, Os and Pt parents, with a view to look for ^°'^Sn and ^^^Sn radioactivities, as well 
as any other new Sn radioactivity with new shell closures in neutrons. Since the cluster 
decays are more probable with daughters as magic nuclei, the decay half-lives are expected 
to drop (be minimum) for the magic daughters. The same idea was utilized earlier for 
the (spherical) sub-shell closed 4oZr daughter [24,25], including also a brief report of the 
results on soSn daughter [24]. This calculation is based on the preformed cluster model 
(PCM) [26,27], described briefly in Sect. II. The results of our calculation are presented 
in Sect. Ill and a summary of our results in Sect. IV. 

2 The model 

The PCM model [26] uses the dynamical collective coordinates of mass (and charge) 

asymmetry, r] = {Ai — A2)/{Ai + A2) and rjz = (^1 — Z2)/{Zi + Z2), first introduced in 

the QMFT [7,8], which are in addition to the usual coordinates of relative separation R 

and deformations ^2% (i = 1, 2) of two fragments. Then, in the standard approximation 

of decoupled R and rj motions, in PCM, the decay constant A or the decay half-life T1/2 

is defined as 

ln2 
A = 7^ = PoPiyo, (1) 

J 1/2 

Here Pq is the cluster (and daughter) preformation probability and P, the barrier 
penetrability, which refer to the rj and R motions, respectively, uq is the barrier assault 
frequency. The Pq are the solutions of the stationary Schrodinger equation in 77, 

h^ d 1 d 
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Figure 1: The decay half-lives Ti/2(s) and other characteristic quantities like preformation 
factors Pq, Q- values (in MeV), and penetrabilities P of different Carbon clusters emitted 
with 5oSn daughters from various isotope of Ba nuclei, calculated on the basis of the PCM, 
plotted as a function of daughter neutron number A'"^!. 



which on proper normahzation gives with a;=0, 1,2,3... . Eq. (2) is solved at a fixed 
R = Ra = Ct{= Ci + C2) (the first turning point of the WKB integral, defined be- 
low), where Ci are the Siissmann central radii Ci = Ri — (l/Ri) (in fm), with the radii 
Ri = 1.28 A,- — 0.76 + 0. 8A^ . Many other radius formulas are available [28] and widely 
used for the calculations of barrier heights, is also a subject of interest for the future study 
in PCM. 

Po = y^,\i^^'>\v{A^))\HVA), (3) 

The fragmentation potential Vii{ri) in (2) is calculated simply as the sum of the Coulomb 
interaction, the nuclear proximity potential [29] and the ground state binding energies of 
two nuclei. 

The proximity potential between two nuclei is defined as 

Vp = 47rC76$(0 (5) 

here 7 is the nuclear surface tension coefficient, C determines the distance between two 
points of the surfaces, evaluated at the point of closest approach and ^{C,) is the univer- 
sal function, since it depends only on the distance between two nuclei, and is given as 
$(^) = -0.5(^ - 2.54)2 - 0.0852(^ - 2.54)^, 

for ^< 1.2511 

(-3.437exp(-^/0.75). /or ^ > 1.2511 
Here, C,= s/b, i.e s in units of b, with the separation distance s=R — Ci — C2. b is the 

diffuseness of the nuclear surface, given by 

7r/2y31n9] (6) 

where tio-^go is the thickness of the surface in which the density profile changes from 90% 
to 10%. The 7 is the specific nuclear surface tension, given by 

2" 



7 = 0.9517 



. N 
1 - 1.7826 



A 



MeVfm-^. (7) 



In recent years, many more microscopic potentials are available that takes care various 
aspects such as overestimation of fusion barrier in original proximity potential, isospin 



effects. A comparison is also available between all models [30]. As noted above, even 
modified proximity potentials were also given. We plan to study cluster decay with these 
new proximity potentials in near future. Here B's are taken from the 2003 experimental 
compilation of Audi and Wapstra [22] and, whenever not available in [22], from the 1995 
calculations of Moller etal. [23]. Thus, full shell effects are contained in our calculations 
that come from the experimental and/or calculated binding energies. We also note that 
for exotic clusters/nuclei with neutron/proton rich matter, new binding energies are also 
available [31]. The momentum dependent potentials and symmtry energy potential which 
are found to have drastic effect at higher densities will not affect decay studies, since 
these happens at lower tale of the density [32,33]. Here in Eq. (4), the Coulomb and 
proximity potentials are for spherical nuclei, and charges Zi and Z2 in (4) are fixed by 
minimizing the potential in r]z coordinate. The mass parameters 5^^(77), representing the 
kinetic energy part in Eq. (2), are the classical hydrodynamical masses of Kroger and 
Scheid [34], used here for simplicity. 

The penetrability P is the WKB tunnelling integral, solved analytically [26] for the 
second turning point Rb defined by V {Rh)=Q,-vaAue for the ground-state decay, and the 
assault frequency uq in (1) is given simply as 

z/o = (2^2//i)'/V^o, (8) 

with E2 = {Ai/A)Q, the kinetic energy of cluster (the lighter fragment), for the Q-value 
shared between the two products as inverse of their masses. Ro is the radius of parent 
nucleus, and /i, the reduced mass. 

3 Results and discussion 

As already stated in the introduction, the cluster decays of various isotopes of seBa to 
78Pt parents are calculated for the daughter nucleus to be always an isotope of 50 Sn nu- 
cleus. For example, for the neutron-deficient ^^*^^^^^Ba and neutron-rich ^^'^"^^"Ba parents 
considered here, different isotopes of Carbon cluster would give rise to various isotopes of 
5oSn daughter. This is illustrated in Fig. 1 for the decay half-life T1/2 of various C-decays, 
together with the Q- values, logarithms of penetrability P and preformation factor Pq, as 
a function of N^,, the neutron number of 5oSn daughter. The impinging frequency Uq is 
nearly constant ~ 10^^ (s~^). All the four quantities Q, P, Pq, and T1/2 show the shell 
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Figure 2: logioTi/2{s) for the most probable clusters emitted from various Ba to Pt parents 
with (a) ^°°Sn (b) ^^^Sn and (c) ^•^^Sn daughter, calculated on the basis of the PCM, plotted 
as a function of cluster proton number Z2. Note the different ordinate-scales used in these 
figures. 



effects at magic A'^d=50 and 82; the Q, P and Pq being large and T1/2 small at these 
numbers. Thus, the most favorable decay is ^^C from ^^^Ba nucleus in the 48 < Nd < 70 
region, leaving behind ^"^Sn as the daughter product, and the ^^C cluster from ^^^Ba in 
the 72 < Nd < 86 region with ^^^Sn as the daughter product. This result is same as 
in Refs. [13] where the most probable clusters for ^°°Sn daughters were obtained as the 
A2 = 4n, N=Z, i^c, 160, 2°Ne, ^^Mg and ^ssi, emitted from the respective Ba to Gd 
parents, and that these were the ^'^C, ^°0, etc., for I'^^Sn daughter, emitted from ^'^^Ba, 
is^Ce, etc. 

In the present study, however, the other most probable clusters considered are (isotopes 
of O, Ne, Mg, Si, S, Ar, Ca, Ti, Cr, Fe and Ni) from heavier neutron-deficient and neutron- 
rich rare-earth parents {'^^-'™Ce, ^'^-^'^Nd, ^^^-^^^Sm, i^s-isogd, i32-i94j3y^ 138-200^1., 
i48-200Yb, i54-208Hf^ 156-208^^^^ i60-2i0Os and i68-2iop^)^ Interestingly, '^C remains to 
be the most favorable cluster-decay from ^^^Ba parent with ^'^'^Sn-daughter [13], but for 
i^^Sn-daughter the most favorable cluster is now ^^Ni from ^^^Pt, instead of ^^C from 
i^^Ba. This is illustrated in Fig. 2(a) and (c), respectively, for ^^'^Sn and ^^^Sn daughters, 
where the most probable clusters emitted from Ba to Pt parents are plotted. The fact 
that the most probable cluster ^*Ni, arising from Pt parents, occur at Nr)=S2 of the soSn 
daughter is illustrated in Fig. 3 for T1/2 alone. However, in Fig. 3, in addition to the 
strong minima at 5oSn-daughter neutrons Nd=82, a new minimum is also shown to be 
persent at No=Q6 for the soSn-daughter, emitting ^^Ni cluster from i^°Pt parent. This 
is further illustrated to be true in Fig. 2(b). Thus, a new possibility of i^^Sn-daughter 
radioactivity is indicated here. Apparently, other cases of interest in Fig. 2 are the ^^Mg 
decay of ^^^Sm and ^°Ca decay of ^^^Yb, respectively, with ^'^^Sn and I'^^Sn daughters. 

Finally, Fig. 4 gives a complete histogram of the decay half-lives logioTi/2{s) as a 
function of the neutron number N^ of the emitted soSn-daughters with the most probable 
clusters (minimum T1/2 values) from some 329 parents taken from Ba to Pt with mass 
numbers A = 110 — 210. Limitation to N^, ~94 since, for N^ >90, the contribution 
from nuclei heavier than Pt would also become important. Note that in Fig. 4, soSn- 
daughter is kept fixed and all possible clusters are considered from different parents (total 
1617 combinations with soSn-daughter and the probable cluster), and then the one with 
minimum half-life time is plotted. Apparently, the shortest half-life time logiQTi/2{s) = 
2.27 (with Q-value=22.16 MeV) is obtained for ^^C decay of ^^^Ba. The role of magic 
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Figure 3: Same as for Fig.l, but for Ti/2(s) alone, and for different Ni clusters emitted 
from various Pt parents. 
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Figure 4: Histogram of logioTi/2{s) versus soSn-daughter neutron number A^^ for the most 
probable clusters emitted from various Ba to Pt parents with soSn as a daughter nucleus 
always, calculated on the basis of the PCM. The associated clusters are also shown on 
the top panel. 
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Nei = 82 is also evident with a minimum in the histogram at ^^^Sn-daughter due to the 
emission of ^^Ni cluster from ^^°Pt parent. The predicted half-life logiQTi/2{s) = 34.974 
(with a Q-value=l 19.292 MeV), which is beyond the limit of present day experiments. 
Thus, as expected, the strongest shell effects occur at N^ = 50 and 82. In addition, 
another minimum due to ^^Ni cluster emitted from ^^°Pt parents could also be of interest 
for a closed shell (either spherical and/or deformed) at No = 66. This minimum is 
comparable to Nd = 82 case, with a predicted decay half-life also of nearly the same 
value {logioTi/2{s) = 34.975, with a Q-value=124.192 MeV), which by all means is again 
very large for experiments. Note from Fig. 2 that the decay half-life for ^^Mg emitted 
from ^^^Sm {logioTi/2{s) = 9.735) lies in between the values for ^^C decay of ^^^Ba and 
^^Ni cluster from ^^*^Pt parent (or ^"^Ni cluster emitted from ^^°Pt parent), rather closer 
to the i^C decay of ^^^Ba. 

4 Summary 

The preformed cluster model (PCM) is used for the cluster decay calculations with 5oSn as 
a daughter nucleus always from various parents of Ba to Pt region. Thus ^"^"^Sn and ^'^^Sn- 
daughter radioactivities is look for the most probable clusters (minimum decay half-life 
time) emitted from the rare-earth parents, and the presence of any new neutron magicity. 
The most probable clusters, respectively, with ^"^"Sn and ^^^Sn daughters, are predicted to 
be ^^C from ^^^Ba and ^®Ni from ^^"Pt parent. Further possibilities with ^"^^Sn and ^^^Sn 
daughters are also noticeable in ^^Mg and ^"^Ca clusters emitted from ^^^Sm and ^^^Yb 
parents, respectively, as the second best new cases. In addition, a new shell is indicated 
at Nj:)=66 with ^^®Sn-daughter due to ^'^Ni cluster emitted from ^^°Pt parents. However, 
at present these calculations seem to be more of an academic interest since the predicted 
half-life times, for at least the ^^^Sn and ^^^Sn-daughter radioactivities, are too large for 
experiments. 

The author is thankful to Prof. R. K. Gupta for many fruitfull discussion. 
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